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A thorough study reveals that the only key parameter for xf [Jlxj), xp') polarization is the ra¬ 
tio (0’^(®S'[®'))/(0’^(®Pg®')), if the velocity scaling rule holds. The ordinary fitting procedure is 
incapable of the determination of this parameter. We provide a universal approach to fixing the 
long-distance matrix elements (LDMEs) for the J/xp and xp' production. Further, with the existing 
data, we implement this approach, and obtain a favorable set of the LDMEs, and manage to recon¬ 
cile the charmonia production and polarization experiment. A quantitative analysis indicates that 
theoretical corrections to the short-distance coefficients change only the values of the LDMEs but 
not the phenomenological results. 
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Introduction. — Nonrelativistic QCD (NRQCD) [l[ is 
one of the most successful effective theory describing 
quarkonium productions and decays (as a review, see 
e.g. i). Despite its contributions, still, it is facing chal¬ 
lenges from many aspects. Three groups [l-@ succeeded 
to accomplish QCD next-to-leading order (NLO) calcu¬ 
lations of J/xp production and polarization at hadron 
colliders, however, with different fitting strategies, ob¬ 
tained quite different values of the long-distance matrix 
elements (LDMEs), consequently leading to different per¬ 
spective of the polarization puzzle. Recently, the LHCb 
Collaboration released their results of rjc hadroproduc- 
tion measurement [^. Three groups [JQ looked into 
the experimental data from relative^ different angles of 
view. Many of the existing works [J, [13, [ill interpreted 
the almost unpolarized experimental results of the J/xp 
hadroproduction measurements as the indication of 
dominance, which violates the heavy quark spin symme¬ 
try regarding the data in Ref. Q. Others SEiil, 
even with different philosophy, also came to the similar 
conclusions. Refs. H, [^ remedied the discrepancy be¬ 
tween the measurements oi J/xp and rjc hadroproduction. 
Notably, their LDMEs are consistent with the velocity 
scaling rule (VSR), which is essential to the NRQCD 
expansion. Even so, both of them failed in the expla¬ 
nation of the J/xp polarization data in midrapidity re¬ 
gions. The polarization parameter A converges in the 
ranges 0.05 < A < 0.2 (for \y\ < 0.6, denoted as El) and 
0 < A < 0.1 (for 0.6 < \y\ < 1.2, denoted as E2) for the 
CMS experiment [l^, and —0.2 < A < 0 for the CDF 
experiment [T^ at y/s = 1.96TeV (denoted as E3). How¬ 
ever, the theoretical predictions [i,iofAforEl,E2and 
E3 reach 0.4, 0.4 and 0.2, respectively. Taking E2 as an 
example, the experimental and theoretical values of the 
ratio of the transverse cross section to the longitudinal 
one (denoted as ^) are about 1.2 and 2.3, respectively. 
A successful effective theory cannot tolerate so large a 
discrepancy. Accordingly, Refs. 0,i both agreed that, 
despite that the yield of rjc and J/xp data were reconciled. 


the corresponding LDMEs still could not solve the J/xp 
polarization puzzle. Interestingly, Refs. @,0 ommit- 
ted the feed down contributions from xp' and Xc to the 
J/xp yield, and claimed that they could ex plai n the J/xp 
polarization data, however. Refs. [E B S llEl indicated 
that, the inclusion of the feed-down parts would ruin the 
results. In sum, the J/xp polarization puzzle is still the 
most challenging question in high energy physics waiting 
for new explorations. The mess of the situation can ac¬ 
tually be attributed to the difficulty in the determination 
of the LDMEs. As is going to be seen later in this paper, 
ordinary fitting procedure is incapable of this question. 
It would be urgent to break through and bring out a 
practical strategy, which is able to definitely either solve 
the J/xp polarization puzzle or phenomelogically disprove 
NRQCD. 

Criticism of the extant fitting strategies. — Before 
putting forward the approach, we first outline the pro¬ 
cedure of determining the LDMEs on the market. Here 
we focus on the direct xp {J/xp, xp') production case, in 
which the cross section can be expressed as [ij da{xp) = 
J2n dfn{0'^{n)), where /„ is the short-distance coefficient 
(SDC) for producing a cc pair with quantum number n, 
and {0'^{n)) is the corresponding LDME. Notice that 
NRQCD is an effective theory, we may expect its pre¬ 
dictions to have an intrinsic deviation (which might not 
be very large, however, does exist) from the reality. In 
addition, our concerns are always limited to specific pro¬ 
cesses (sometimes because of the lack of knowledge on 
other processes, which is due to e.g. experiment is lack¬ 
ing or higher-order corrections are large). For this reason, 
we regard two sets of the LDMEs leading to close pre¬ 
dictions in the processes we concern as ’’equivalent for 
these processes”. Further progress in both theoretical 
calculation and experimental measurement would distin¬ 
guish the ’’equivalent” sets of the LDMEs. Up to QCD 
NLO, perhaps the xp hadroproduction is the only process 
in which the dominant contributions are all counted. 

When we fit the xp yield data, the standard deviation 
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is a quadratic function of the LDMEs. Here, , cr“ and 
ed denote the theoretical prediction, and the experimen¬ 
tal central value and error for the dth experimental data 
point, respetively, and D is the degree of freedom in the 
fit. By way of illustration, we only take the three color- 
octet (CO) matrix elements, 

and (0’^(3p^®')), as to be determined. To keep the homo¬ 
geneity of the dimensions of the CO LDMEs, in this pa¬ 
per, we definde /jpis] and by multiplying and 

dividing by a factor of respectively. For convenience, 
(O’^(n)) is alternatively abbreviated to x 10“^ GeV^ 
in the following, with n =1,2,3 representing 
and respectively. 

_Ordinary fitting procedure is to solve the equations, 
dx^jdOn = 0, to fix the values of the LDMEs at which 
the reaches its minimum. However, Ref. [13 found 
that the SDCs for the three CO channels roughly satisfy 
a linear relation 


^o/*ic[8] “b J^[8] , (2) 

O ‘^0 ‘^1 

thus, only two of the three LDMEs can be fixed through 
the fit of the yield data. For instance, the cross section 
for direct ip hadroproduction can be expressed as 

daW = Agis] Mf, (3) 

where and Mf are defined by 

Mf = {O'^Csf^))+ri{0'^(^pf^^)). (4) 

One can fit the yield data and obtain the values of Mq 
and Ml by employing Eq(|31). 

The reduction strategy provided in Eq.Q is feasible to 
work on the ip yield, nevertheless, we find that it is not 
suitable for the polarization problem. On the one hand, 
A is sensitive to as many as two parameters, namely Mi 
and ri; even a slight variation of the two parameters 
can cause dramatic change of A. On the other hand, 
Eq.(I5]) is only an approximate relation; vq and ri are 
different in different kinematic regions. For instance, for 
the CDF experimental condition, in the range 7 GeV < 
Pt < 30 GeV, To = 3.9 and ri = —0.56, while in the 
range 11 GeV < p* < 30 GeV, rp = 3.5 and ri = —0.53. 
The difference of ri for the two pt ranges is large enough 
to completely change the predictions of the polarization. 

New methodology and its implementation .— Ref. @ 
provided an evidence for the VSR, which is the most 
fundatmental bases of NRQCD (otherwise, the infinite 
higher excited Fock states of cc will be involved). To 


this end, we constrain our discussions within the ex¬ 
tent where this rule is not violated. Under this assump¬ 
tion, channel cannot dominate the ip production, 

thus, the and ^Pf^ channels must contribute a large 
part and the cancellation between them would be signif¬ 
icant. We emphasize that, the cancellation is actually 
natural, since, at QCD NLO or higher order, the two 
channels are associated; only the combination of the two 
channels is divergence free and NRQGD scale indepen¬ 
dent. We need only to argue that for any process, the 
cancellation does not cause unphysical results (negative 
cross sections). Actually, it is a direct conclusion of the 
fragmentation mechanism. The combination (as well as 
the cancellation) of the two channels can be carried out 
at fragmentation-function level. In addition, since the 
channels cannot saturate the ip hadroproduction, 
the leading and next-to-leading power [isl . Il^ terms in 
the combination of the and ^Pf^ channels must be 
positive. Gonsequently, the combination of the fragmen¬ 
tation functions (FFs) multiplied by the LDMEs of the 
two channels must be positive definite. So what we need 
to work out is how much do they cancel, or equivalently, 
what is the value of 

We will find that the ip polarization is extremely sen¬ 
sitive to while the ip yield is not. This brings us to a 
subtle circumstance that small variations of result in 
equivalent LDMEs for the ip production, however, give 
totally different predictions of the polarizations. This is 
the exact reason why one cannot succeed in the explana¬ 
tion of the Ip polarization by employing the LDMEs ob¬ 
tained in the fit of the yield data by minimizing the x^- 
Since for the J/ip, the LDME obtained in Ref. Q 
has a large uncertainty, while for the ip', it is totally un¬ 
known, in our strategy, we assign it different values, and 
directly fit (0’^(^S'[®^)) and (0’^(^Pg®^)), and find that the 
variation of (0 '^(^S'q®^)) only leads to equivalent LDMEs 
for both Ip yield and polarization. Having this, we can 
assign it any possible value under the constraint of the 
VSR and Ref. Q. 

In low and high pt region and forward (backward) ra¬ 
pidity (y) region, large logs {log{ml/pf), log{pt/E^)) are 
important. We constrain our concerns in the kinematic 
region, \y\ < 1.6 and 7 GeV < p* < 30 GeV for J/ip 
and 11 GeV < pt < 30 GeV for ip' [I^, in order to 
keep the perturbative expansions safe. For the yield, all 
the data in these region provided in Refs. [20l-l^ are in¬ 
cluded in our fit. Note that the ip' polarization data in 
Refs. [Il[il have large errors, we would start from the 
study of the J/ip, where, for the yield, the contributions 
from Ip' feed down are counted, while for the polarization, 
they are neglected, under the consideration that ip' feed 
down only contribute a fraction of less than 10% to the 
prompt J/ip production, which cannot affect the J/ip po¬ 
larization. After acquiring some educational knowledge 
from the J/ip case, we would come back to study the ip' 
meson. 

For the Xc production, we adopt the same parameter 



3 


choices as in Ref. [1^. To calculate its contributions 
to the polarized J/ijj) we follow the scheme developed 
in Ref. The parameters adopted in the analysis of 
J/tp and Ip' yield and polarization are identical to those 
adopted and obtained in Ref. @ , except for the values of 
(qV-( 35'[81)) and Since the color-singlet (CS) 

LDME and are identical to those in Ref. 

we can expect that they are naturally consistent with the 
rjc hadroproduction data Q. 


CO 

o 

1.50 

1.60 

1.70 

1.80 

1.90 

2.00 


0.898 

0.934 

0.971 

1.008 

1.044 

1.081 

Rj/t/j 

0.599 

0.584 

0.571 

0.56 

0.549 

0.540 


2.16 

2.03 

1.98 

2.00 

2.10 

2.27 


TABLE I: The value of and the corresponding Rj/^ 

and at each specific value of . The global error of 
is ±0.011. 


We hrst directly fit the ip' yield and obtain Of = 
0.48 ± 0.02 and of = 0.80 ± 0.05, where of = 0 
is set as default. Employing them, associted with the 
Xc predictions ,we can extract the direct part from the 
prompt J/tp yield data, and directly fit the data and ob¬ 
tain = 1.0 ± 0.1 and = 1.7 ± 0.1, which are 
consistent with our previous work Q, with ra 1.98. 
We remember that small deviations from the optimized 
values of the LDMEs provide equally good descriptions 
of the yield data, for this reason, we fit at each 

specific value of . The results are listed in TABlH 
Regarding Eq. o, the deviation of the yield curve for 
= 2.0 from that for = 1.7 is less than 10% 
of the experimental error in average. Higher order cor¬ 
rections, large log resummations, experimental errors, or 
even numerical uncertainties and the intrinsic errors of 
an effective theory are comparable with that. In other 
words, the LDMEs listed in TABU] are equivalent for the 
J/ip yield. We can summarize the LDMEs in TABjTjin a 
compact form as 


Ot = k^Ot + (5) 


while for J/ip^ 

kj/^ = 0.367, bj/^ = 0.348 ± 0.011. (6) 

Employing Eq.@, we can fit the J/ip polarization 
data. The degree of freedom has been reduced to one. 
Since the measurement of CDF Run 1 and Run 11 are 
contradict with each other, we give up using the Run 1 
data for their large uncertainties. Considering that, the 
polarization is a ratio; even a slight error (as small as 
20%) can cause significant deviation, we drop the data 
in low pt region, where the precision provided by the 
perturbative expansion is quite difficult to control. Only 
the Pt > 10 GeV data in Refs. 013 are adopted in our 




FIG. 1: J/ip and ip' yield at the Tevatron and the LHC in 
medium Pt region. All the LDMEs follow Eq.(|5]) in association 
with Eq. for J/ip and Eq.Q for ip'. The data are taken 
from Refs. [20]4^ 


fit. Including the contributions from Xc feed down and 
excluding those from ip' ^ we obtain the value of as 

= 0.546 ±0.006. (7) 

We emphasize again that Rj/^ is the only parameter to 
govern the J/ip polarization, as long as the VSR is kept. 
For instance, if we fix Rj/^ and vary from 

its upper to lower bound obtained in Ref. Q, or vary 
from 1.5 to 2.0, the corresponding change of A is 
less than 0.02 (most of the time, much smaller than this). 
Accordingly, Eq.® and Eq.® provide the uncorrelated 
form of the LDMEs; the uncertainties of and Rj/^ 
describe those of the J/ip yield and polarization, respec¬ 
tively. 

The L.H.S. plot of Fig |T] displays the theoretical pre¬ 
dictions versus data for the J/ip hadroproduction. The 
bands are expanded by the curves for 1.7 < Og^’^ < 2.0, 
with the corresponding obtained through Eq.® 

and Eq.®. The LDMEs in this range can provide 
equally good descriptions of the yield data. However, as 
is shown in Fig® they result in totally different polar¬ 
ization predictions. The solid curves are produced with 
the LDMEs obtained in Refs. [1,0, which corresponds 
to Og^’^ = 1.7 with Eq.® and Eq.® satisfied, while 
the bands are produced with the LDME ranges obtained 
in Eq.®, which corresponds to 1.88 < Og^’^ < 2.01. 
We can see that the CMS data are well described in our 
framework. For the CDF data, the discrepancy between 
theory and experiment is larger. We need to present 
the values of ^ to acquire a better understanding of the 
problem. For the pt = 20 GeV data, ^ = 1.2 and 0.7 for 
the theoretical prediction and the experimental central 
value, respectively. A correction of about only 20% can 
fill this gap. In the next part, we will further argue that, 
when Pt < 20 GeV, the polarization is sensitive to the 
corrections while at higher pt, it is not. In this sense, we 
can say that, the J/ip yield and polarization data as well 
as the Pc yield data are reconciled within the NRQCD 
framework. 

Note that the CMS data for the ip' polarization is not 
monotonic with respect to the rapidity [l^ . we can ex¬ 
pect that precise measurement will significantly change 
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FIG. 2: J/ijj polarization at the Tevatron and the LHC. 
The solid curves are produced with the LDMEs obtained in 
Ref. [^, while the bands are produced with the LDMEs cor¬ 
responding to Eq. 0 . The data are taken from Refs. MM- 


Ot 

0.50 

0.60 

0.70 

0.80 

0.90 

1.00 

1.10 

of 

0.353 

0.395 

0.437 

0.479 

0.521 

0.563 

0.605 

R^i)' 

0.706 

0.658 

0.624 

0.599 

0.579 

0.563 

0.55 


1.23 

0.79 

0.53 

0.44 

0.53 

0.79 

1.22 


TABLE 11: The value of O 2 and the corresponding and 
at each specific value of Of , while Of = 0 is set as 
default. The global error of O 2 is ±0.003. 


the central values. For this reason, we would see how 
much our prediction of the polarization can reach, while 
the yield data is well described, instead of carrying out 
a fit of the experiment. We can obtain a similar table, 
TABini as TABUfor J/ip, and the corresponding linear 
relation with 

V = 0.42, 5^,' = 0.143 ± 0.003. (8) 

The bands in the R.H.S. plot of FiglT]and in Figl3] cor¬ 
respond to the range 0.6 < Of < 1 . 0 , in which, as is 
displayed in TABUTI the s are quite small. This range 
of the LDMEs lead to quite large polarization bands, 
yet, not large enough to cover all the data points, as is 
shown in FigEl Therefore, we also present the curves 
for Og = 1 . 1 , which can cover the upper bound of 
the error bands of the CMS data in the rapidity range 
0.6 < |y| < 1.2. The curves for the ip' yield also are in 
good agreement with the data, albeit the x^’s are rela¬ 
tively larger. Figl3] and Fig0 clearly manifest the fact 
that, when Eq.Q is held, the yield data can be well re¬ 
produced in very large range of the LDMEs, while the 
polarization is extremely sensitive to . This is to say, 
it is almost impossible to describe the polarization data 
using the LDMEs obtained through the fit of the yield 


FIG. 3: ijj' polarization at the Tevatron and the LHG. All the 
LDMEs follow Eq.0 and Eq.0. The data are taken from 
Refs. P,[T^. 


data; even the variation of the yield curve is as slight as 
the intrinsic error of an effective theory, the polarization 
will change dramatically. We can also conclude that, it 
is almost certain that, when higher-precision polarization 
data come out, the yield and polarization of tjj' can both 
be well reproduced. 

Having got O 2 and Og for Of = 0, we can attempt 
to assign Of a larger value consistent with the VSR. 
When of = 1.0, we obtain k^i = 0.42 and = 0.115± 
0.002. Following the same procedure, we find that, the 
phenomenological results does not change, which proves 
that varying Of and redoing the fitting procedure only 
leads to equivalent LDMEs. 

Universality .— Since the and channels suf¬ 

fer from large cancellations, we need to further discuss 
the stability of the results, or equivalently, to answer the 
question, whether a slight (or even large) correction can 
ruin the phenomenology, or equivalently, whether the so¬ 
lution of the LDMEs keeping the polarization lies in the 
region where the for the fit of the "0 yield is not too 
large. 

Most of the time, the corrections are proportional to 
the NLO results. Without lost of generality, we assume 
that the unpolarized SDCs remain, while only the ratio 
of the transverse part to the longitudinal one and 

Cspisi) changes. Under this assuption, the linear rela¬ 
tions in Eq.® hold, so do TABU and TABITTI Only the 
relation between A and changes. We recalculate the 
values of A for specific values of ■Cagisi and ^ 3 p[ 8 ] as a func¬ 
tion of R^ numerically, and find that, in large pt region, 
especially pt > 20 GeV, to obtain the previous results of 
the polarizations, the corresponding x^’s for the fit of the 
ijj yield are also small. For instance, if decreases by 

30% (which is quite a large correction), to keep the po- 
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larization of the J/ip or ip' invariant, need to increase 
0.05 at the most. For both J/ip and ip', the for the fit 
of the yield data at the shifted is tolerable regarding 
TABUand TABITTI However, in the region pt < 20 GeV, 
A is more sensitive to the corrections. Most of the time, 
the corrections would enhance the longitudinal fraction, 
so, we can expect that, the CDF data for the J/ip polar¬ 
ization will be better described when further corrections 
are counted. 

We need also to investigate whether high pt data can 
help to fix the LDMEs in the absence of the polarization 
data, if the FFs dill are employed. Actually, both 
the and SDCs scale as pp'^ in large pt limit, 
thereafter, they will always be tangled. For instance, the 
ratio of the curve for Of = 0.395 and Of = 0.60 to 

the curve for = 0.605 and O 3 = 1.10 is about 1.5 
in large pt limit in midrapidity region, while the differ¬ 
ence of the corresponding A between the two sets of the 
LDMEs can be as large as 0.7. Therefore, high pt data 
might provide slight constraints of the LDMEs, still, it is 
impossible to describe the polarization with the LDMEs 
obtained employing only the yield data. 

The only remaining question we would discuss is that, 
whether a global fit, employing both the yield and po¬ 
larization data, can give the same LDMEs as presented 
in Eq.® and Eq.([7|). First we set the LDME for 
as a free parameter, and obtain = 3.0 ± 1.6, 

02^^ = 0.7 ± 0.2, and = 1.3 ± 0.4, which con¬ 


tradicts with the r]c hadroproduction data. So, we set 
= 0.78 according with Ref. [^. The global fit gives 
= 1.02 ± 0.03 and £> 3 /’^ = 1.84 ± 0.06, the corre¬ 
sponding value of Rj/^jj of which is Rj/ip = 0.554, which 
is above the upper bound of the range in Eq. ©• So, we 
can conclude that, a global fit, even including the po¬ 
larization data, still cannot tackle the J/ip polarization 
puzzle. 

Summary .— In this paper, we discovered the unique 
key parameter which governs the ip polarization, namely 
RtP = and reconciled all the char- 

monia yield and polarization data in midrapidity region 
within NRQCD framework. When i ?.0 is fixed, varying 
the LDMEs hardly changes the polarization if the VSR 
is not violated. Besides, we found that the polarization 
is extremely sensitive to even under the constraint 
of the yield data. Accordingly, it is almost impossible 
to explain the polarization with the LDMEs fixed by the 
yield data. Through a brief analysis, we found that the 
cancellation between and is natural, and fur¬ 
ther corrections might change the values of the LDMEs 
but not the phenomenological results. A global fit, even 
including the polarization data, is incapable of the deter¬ 
mination of the LDMEs. 

We thank Yan-Qing Ma for helpful discussions. This 
work is supported by the National Natural Science Foun¬ 
dation of China (Nos. 11405268). 
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